function S=area(t,T)
%积分计算阴影部分面积
    [~,k1]=max(T);
    k=find(T>=217);
    k0=k(1);
    dt=diff(t);
    dt=dt(1);
    DT=T-217;
    S=sum((DT(k0:k1-1)+DT(k0+1:k1))/2*dt);
end

